Photo-induced insulator-metal transition of a spin-electron coupled system 
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The photo-induced metal-insulator transition is studied by the numerical simulation of real-time 
quantum dynamics of a double-exchange model. The spatial and temporal evolutions of the system 
during the transition have been revealed including (i) the threshold behavior with respect to the 
intensity and energy of light, (ii) multiplication of particle-hole (p-h) pairs by a p-h pair of high 
energy, and (iii) the space-time pattern formation such as (a) the stripe controlled by the polar- 
ization of light, (b) coexistence of metallic and insulating domains, and (c) dynamical spontaneous 
symmetry-breaking associated with the spin spiral formation imposed by the conservation of total 
spin for small energy-dissipation rates. 
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Changes of electronic states induced by the photo- 
irradiation have been the subject of intensive studies 
in chemistry, biology and physics. For example, the 
photo-induced chemical reactions in molecules are the 
basis of many biological functions. In solids, the struc- 
tural changes induced by the absorption of photons and 
the relaxation after it are a common phenomenon [l|. 
In particular, in correlated electrons such as the tran- 
sition metal oxides, the collective nature of the system 
enhances the sensitivity to the external stimuli includ- 
ing photons, and hence even the weak photo-irradiation 
can trigger the changes of whole system as the phase 
transition 0-ly|- This is because in strongly correlated 
electronic systems, the rich phases of spin/charge/orbital 
orderings with nearly degenerate energies compete with 
each other [7|. For example, in the case of mangan- 
ites, colossal magneto-resistance effect occurs because the 
weak external magnetic field drives a transition from the 
insulating charge/spin/orbital ordered state to the ferro- 
magnetic metallic state [8| , and similar transition can be 
triggered also by the photo-irradiation Q. 

In contrast to these ample examples of photo-induced 
phase transitions, the theoretical studies arc still in the 
premature stage. Although a classical model for the 
photo-induced phase transition has been proposed [9|, 
the quantum theory of photo-induced dynamics of corre- 
lated electrons are limited to small size systems or to one- 
dimensional systems p^l - [l2| . Recently, we have proposed 
a model for the real-time quantum dynamics of corre- 
lated electrons, i.e., a double exchange model [13 where 
the electrons are interacting with the classical spins, to 
simulate the relaxation after the photo-excitation |14| . 
In that work, wc have developed a theoretical method 
for the fully quantum-mechanical time-evolution of elec- 
tronic wavefunction combined with the Landau-Lifshitz- 
Gilbert (LLG) equation for classical spins, which enables 
the analysis of the larger-scale/higher-dimensional sys- 
tems. By that theoretical study, the self-organized space- 
time structure induced by the quantum transitions dur- 
ing the relaxation has been revealed |14| . 



In this paper, we present the real-time quantum dy- 
namics of the photo-induced insulator-to-mctal (IM) 
transition for a generalized double exchange model show- 
ing the competing antiferromagnctic insulating (AFI) 
and ferromagnetic metallic (FM) states separated by the 
first-order phase boundary. Many unexplored issues can 
be addressed in this model, i.e, the threshold behavior 
with respect to the intensity/energy of light, the effect 
of the polarization of light, space-time pattern formation 
during the transition, the role of the conservation of total 
spins, and the asymmetry between the two directions of 
the transition. 

Wc start with the Hamiltonian, 



H ^ -t '^ cl^Cjs + h.C. - Jh"^ 4s^is'^ss' ■ Si 
<ij>,s iss' 



<ij> 



<ij> 



where <ij> denotes a nearest-neighbor pair, s and s' 
are indices for electron spin, respectively, and t? = 
(cr^, cr^, cr^) are Paidi matrices. The local spins, S'i's, 
are taken to be classical vectors with magnitude S, and 
other notations are standard. We consider the half-filled 
case, i.e., one electron per site. In Eq.([T|), we have intro- 
duced the biquadratic exchange interaction with coupling 




FIG. 1: (a) Ground state phase diagram of the model Eq. ([T| 
in the plane of S^J and SJh with fixed 5* Jjv=— 0.2 in units of 
t. There exist three phases, i.e., ferromagnetic metal (FM), 
antiferromagnetic insulator (AFI) and canted antiferromag- 
netic insulator (CAFI). (b) Lowest energy per site measured 
from ground state as a function of angle between neighboring 
local-spins in a two sub-lattice consideration on 8 x 8 system 
(see text). 



constant Jn, which gives an energy barrier between fer- 
romagnetic (F) and antiferromagnetic (AF) states when 
Jn is negative. 

Figure 1(a) shows the ground state phase diagram of 
the model Eq. ([1]) on the two-dimensional square lattice, 
in the plane of S^J and SJh with fixed S'^Jn=— 0.2 in 
units of t. The solid curve indicates the boundary of 
first-order phase transition. In the right side, the anti- 
ferromagnetic spins open a gap and the system is insu- 
lating. In the left side, the spins are perfectly or almost 
ferromagnetic. In the ferromagnetic state, the system 
is metallic, but the small staggered component opens a 
small gap for SJh^^-2 ^M,- The phase transition at 
SJh— ^-2 (broken line) is of the continuous second order. 

Let us focus on the boundary of first-order phase tran- 
sition. Figure [Ijb) shows the energy as a function of the 
angle 6 between the sublattice spin directions for 8x8 
size system in periodic boundary condition with a pa- 
rameter set, t=l, SJh=-1, S^J=-0M3, S'VAr=-0.017, 
S—1. This condition is used hereafter. The AF state and 
the F state, respectively given by 6' = tt and 0, are al- 
most degenerate, and show the local stabilities with the 
potential barrier between them. 

We introduce the Landau-Lifschitz-Gilbert (LLG) 

equation for the motion of the local spins. Si = /ie//.j x 

S, - aS, X S,, where /J^//,, = (*(r)| - dH/dS,MT)) 

with \^{T)) being the electronic wavefunction of a 
Slater determinant at time T obtained by solving the 
Schrodinger equation as discussed in Ref. [ij]. The 
Gilbert damping constant Sa describes all the other re- 
laxation processes of the spins than the coupling to the 
conduction electrons. 

Figure [5] shows time-evolution of the photo-excited 
electronic state and its relaxation dynamics for Sa ~ 1. 
The unit of time T is h/t and hence is typically ~ 10^^^ 
sec assuming t = 0.4eV. We start with the AFI state, i.e., 
the electronic state has an insulating energy-gap 2SJh 
and lower energy band is occupied by electrons and the 
upper energy band is empty at T = 0. In order to mimic 
the thermal fluctuation, we introduce a random tilting 
of each spin from the perfect AF configuration up to 0.1 
rad which corresponds to the state with an excitation 
energy of ^ O.OOOSi from the ground state. The photo- 
excitation is achieved by introducing the time-dependent 
vector potential into the hopping matrix element: for 
the nearest-neighbor pair <ij> in y-direction, we use 
t — > te*"^^-^' where the phase A{T) of the hopping ma- 
trix element is given by the vector potential. In the early 
stage, T = - 80, we have applied A{T) = O.lsin(wr) 
in this simulation. The frequency w is tuned for the en- 
ergy difference between the second highest and the sec- 
ond lowest energy levels of the AF state on the finite size 
system. The particles are excited from lower to upper 
energy band (see Fig. [2ja)) and the particle- hole (p-h) 
pairs are created on the second highest and the second 
lowest energy levels (see Fig. [^Jb)). In the excited state, 
the spatial distribution of excitation energy is uniform as 
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FIG. 2: (Color online) Time evolution of the photo-excited 
electronic state and its relaxation dynamics for Sa =1. (a) 
Black line with scale on left axis is the time dependence of to- 
tal occupation number above Fermi level, riex- Red (gray 
in print) line with scale on right axis shows the time de- 
pendence of the equal-time spin-structure-factor, defined by 

S{q) = ■^Y.,j^^ ■ Sje^'^i'-'), for q = (7r,7r), where i and 
j are the vector representation of the lattice points i and 
j, respectively, q is that of the reciprocal lattice-point, and 
A^ denotes the system size. The arrows on top axis are the 
marks for T = 80, 152, 430 and 1000, respectively, (b) Elec- 
tron occupancy at T = 80 represented by the black area with 
the density of states (DOS) (thin line). The Fermi energy 
is taken to be zero. In the same way, the electron occu- 
pancy at T = 152, 430 and 1000 are shown in (c), (d) and 
(e), respectively. Figures (f), (g) and (h) are the distribution 
of excitation-energy-density in real space measured from the 
ground state at T — 80, 152, and 430, respectively. Numbers 
on horizontal (vertical) axis is for the coordinates of lattice 
points in x (y) direction, (i) S{q) at T = 80. (j) Transverse 
component of the equal-time spin-structure-factor S±{q) at 
T = 152. (k) The same with (i) but T = 430. A spline 
interpolation is used for figures (f)-(k). 

seen in Figl2][f). At T = 80, total amount of the excited 
electron is about 4, and after that, A{T) = 0. 

Up to r ~ 200, the local-spin structure is almost AF 
one, so that the density of states (DOS) hardly changes 
as shown by thin lines in Fig.[2Ub) and (c). However, the 
electronic state has been changed dramatically, i.e., the 
electron occupation shows a time evolution and the spa- 
tial distribution of excitation energy forms a stripe pat- 
tern as seen in Figs. [2Jc) and (g). Correspondingly, the 
transverse component of the local-spins shows a dynam- 
ical pattern as in Fig. [2IJ). During the photo-excitation 
by the light polarized in y direction, the electrons are 
accelerated along the direction. The polarized motion of 
electrons derives a deformation of the local-spin struc- 
ture to be favorable for the electron motion. The charac- 



teristics of the relaxation dynamics appear in the stripe 
pattern as seen in Fig. [H^g) reflecting the wave function 
and the spatial distribution of the excited p-h pairs. 

In the period T = 200 ^ 300, the system shows 
a drastic change, that is, AFI to FM transition. For 
T < 200, the spin-structure-factor S'(7r,7r) is almost 1, 
and we clearly see the insulating energy-gap 2SJh be- 
tween upper and lower energy bands, i.e., the system 
is the (excited) AFI. For T > 300, on the other hand, 
S'(7r,7r) is almost 0, and the gap closes. In fact, the en- 
ergy level structure is very close to that of the tight- 
binding Hamiltonian for free-electrons with a Zeeman 
splitting. (The sparse energy level distribution around 
zero in energy axis is due to the finite-size effect.) As 
seen in Fig. [2|k), the local-spin becomes the F state cor- 
responding to the metallic electronic-states. The elec- 
trons, however, remain still highly excited. As seen in 
Figs.[2](a) and (d), the total number of excited electrons 
above Fermi level increases and shows a peak at T '^ 430. 
After that, it decreases but is as large as the number of 
initially excited electrons, up to T ^ 660. On the other 
hand, it has been confirmed that the total energy of the 
system, i.e., the electrons and the local spins, decreases 
monotonically due to the Gilbert damping. 

The enhancement of the total number of excited elec- 
trons above Fermi level is due to the Auger process. This 
is a characteristics of the interacting electron systems, 
and is in sharp contrast to the relaxation dynamics of 
conventional semiconductors which physics is well de- 
scribed by the single carrier problem. This multiplica- 
tion of p-h pairs during the relaxation dynamics offers 
an interesting possibility to enhance the photo-current 
generation and solar cell reaction. 

This AFI-to-FM transition, however, does not occur 
when the excitation is not strong enough. For the results 
shown in Fig. [2l the photo-excitation has been done in 
the early stage, < T < T/ = 80, with w ^ 7SJh- By 
decreasing the period T/, we can reduce the number of 
the p-h pairs by the photo-excitation nex{Tf). We find 
that the IM transition hardly occurs for Ue^iTf) < 2, 
in the present numerical condition. 

We have also examined the frequency w dependence 
for the IM transition. In the case that w = 2SJh, the 
threshold of the IM transition is lying on 4 < nex{Tf) < 5 
in the numerical simulations. Although the IM transition 
occurs even in this case, we do not find the Auger pro- 
cess. Comparison with the above case of w = 7SJh, it is 
concluded that the higher energy photon is more effective 
to induce the IM transition through the Auger process. 

When nex{Tf) is close to the threshold, we find the 
coexistence of the insulating and metallic domains. The 
results for u! = 2SJh and nex{Tf) ~ 4.13 are summarized 
in Fig. [31 The local-spin state changes gradually in a long 
time period, 150 < T < 1100 (see the time dependence 
of ^(Tr, tt) in Fig.IS^a)). A long time later, the AF (blue) 
and F (yellow) domains arc spatially separated as seen 
in Fig. [3] (d), and the corresponding excitation energy 
density is shown in Figs. |3] (e). The energy levels and 
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FIG. 3: (Color online) (a) Black line with scale on left axis is 
the time dependence of total energy per site measured from 
ground state. Red (gray in print) line with scale on right axis 
shows S{n, n) as a function of T. For the photo-excitation, 
A{T) = 0.15sin(2SJ//T) is used. Different time scales are 
used for T < 5 and 5 <T, respectively, to zoom up the initial 
stage for the photo-excitation. During the photo-excitation, 
total energy increases, and later on (5 < T) it decreases due to 
the Gilbert damping, (b) The same with Fig. [g^b) but T = 6 
in this simulation, (c) The same with (b) but T = 2000. (d) 
Local spin structure at T = 2000. The stick with dot indicates 
the local spins. There exist AF and F domains. The F and 
AF domains are expressed by yellow (gray) and blue (dark 
gray) regions, respectively. In the figure, the spin belongs to 
F (AF) domain when all the neighboring spins are in F (AF) 
alignment, and the other spins are put into the white area. 
(e) The same with Fig. ^i) but T = 2000 in this simulation. 

occupation numbers before and after the transition are 
shown in Figs, ^h) and (c), respectively. The residual 
p-h pairs in Fig. [3] (c) are confined in the F metallic 
domain, so that the electronic state form a slightly higher 
energy density region on the F domain and the stable AFI 
domain (see Figs. [3l^d) and (e)). This meta-stable state 
continues for a long time within our simulation (at least 
up to T - 6000). 

We have also studied the more realistic case of 
5q;=0.01 (Fig. 21). Other parameters are the same with 
the case shown in Fig. [51 Due to the reduced relax- 
ation rate, it takes much longer time for whole relax- 
ation dynamics compared with that in Fig. [51 Despite 
the reduced relaxation rate, the real-time dynamics of 
the electronic and local-spin states is not silent and is 
rather active. During the relaxation dynamics, the direc- 
tion of the staggered magnetization changes drastically. 
And the spatial inhomogeneity of the electronic state has 
also been developed (sec Fig. [^f)). The slow energy re- 
laxation is advantageous for the Auger process resulting 
in the large amount of p-h pair creation at T ^ 8000 
when the gap closing occurs (see Fig.[31Ja)). At this tran- 
sition, all the local-spins play very active dynamics with 
developing spatial inhomogeneity of the electronic state 
(see Fig. [H^g)). The change of spin configuration from 
AF state is clearly seen in the equal-time spin-structure- 
factor shown in Figs.[31Ja) and (j): The intensity at (tt, tt) 
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FIG. 4: (Color online) The same with Fig. [2] but Sa = 0.01 
is used. In (a), different time scales are used for T < 80 
and 80 < T, respectively, to zoom up the initial stage. The 
arrows on top axis are the marks for T — 1000, 8110 and 
20000, respectively. Figure (b) is the same with Fig. [2jb) 
but Sa — 0.01. The figures (c), (f) and (i) show the results 
at T=1000. The figures (d), (g) and (j) show the results at 
T=8110. Here, S{q) is presented on (j), instead oi S±{q). The 
figures (e), (h) and (k) show the results at r=20000. 

becomes small and that at around the center region is 
grown. Near the center, there exist four spots, and two 
of them disappear later on in the time evolution. Fi- 
nally, a spiral spin-structure is obtained in this case (see 
Fig. Hl^k)). Due to the Gilbert damping Sa, the total 
spin-angular-momentum can change in the relaxation dy- 
namics. For large Sa the spin state can be easily changed 
from AF to F as in Fig. ^g) but not for small Sa be- 
cause the conservation law becomes effective. Therefore, 
the ferromagnetic state is avoided and the spiral spin con- 
figuration appears characterized by the equal-time spin- 
structure- factor shown in Fig. |U[k) . In the present case, 
the propagation vector of the spiral spin-state is in x di- 



rection. This is the dynamical symmetry breaking, and 
the direction {x in this particular case) is chosen by the 
accidental reason while the memory of the polarization 
of photo-excitation {y direction) has been already lost. 

So far, we have discussed the AFI-to-FM transition. 
As shown in Fig.[TJb), FM and AFI states are almost de- 
generate, so that FM-to-AFI transition of this system is 
worth to be examined. In the F state, however, the exci- 
tation with g = is prohibited. Therefore, we have pre- 
pared the initial excited-state within F state by putting 
the initial electron occupation by " hand" . Although we 
have examined a number of numerical simulations, mctal- 
to-insulator transition has not been obtained. The nu- 
merical results suggest that metal-to-insulator transition 
is much more difficult than insulator-to-metal transition 
in the present model. 

In summary, we have studied the insulator-to-metal 
transition by the photo-irradiation in a model of com- 
peting AF and F states, (i) The threshold behavior with 
respect to the intensity and energy of lights, (ii) Auger 
process, i.e., multiplication of particle- hole (p-h) pairs 
by a p-h pair of high energy, (iii) the space-time pattern 
formation including the stripe controlled by the polariza- 
tion of light, and (iv) the intriguing spin spiral state by 
the conservation of total spin for weak Gilbert damping 
have been revealed. The photo-carrier generation under 
the external electric field and its implications to the solar 
cell reaction arc the issues left for future studies. 
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